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Abstract: Recently a remarkable agreement was found between lattice simulations 
of long Wilson lines and behavior of the Nambu Goto string in fiat space-time. 
However, the latter fails to fit the short distance behavior since it admits a tachyonic 
mode for a string shorter than a critical length. In this paper we examine the question 
of whether a classical holographic Wilson line can reproduce the lattice results for 
Wilson lines of any length. We determine the condition on the the gravitational 
background to admit a Coulombic potential at short distances. We analyze the 
system using three different renormalization schemes. We perform an explicit best fit 
comparison of the lattice results with the holographic models based on near extremal 
D3 and D4 branes, non-critical near extremal AdS^ model and the Klebanov Strassler 
model. We find that all the holographic models examined admit after renormalization 
a constant term in the potential. We argue that the curves of the lattice simulation 
also have such a constant term and we discuss its physical interpretation. 
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1 Introduction 



Recently lattice simulations of Wilson lines were found to admit a remarkable agree- 
ment with the relation between the energy and the length of the string as follows 
from the Nambu Goto (N.G.) formulation. The simulations measured this relation 
in two [1] and three [2] space dimensions for both closed strings wounded along a 
circle as well as for the open string case [3]. The dependence of the energy on the 
distance between the endpoints of an open bosonic quantum string in d dimensional 
fiat space-time is given by [4] 



E=\j{a^Y + {n-^)^1ra (1.1) 

where a is the string tension, I is the separation distance and n is the excitation 
level. It was found out that for the ground state the agreement with the full square- 
root formula is better than the agreement with any approximation of this formula 
for long strings {af » 1), like E ~ a£ + + + -• • What is in 

particular remarkable is the fact that the agreement holds also for short distances, 
up to ly/a ~ 1.6. 

Note however that the N.G. behavior (1.1) cannot serve as an adequate descrip- 
tion for distances shorter than the critical length 

^ ' 1.2 



since there the energy develops an imaginary part. It is thus clear that to accom- 
modate the short range behavior of the potential one needs to go beyond the stringy 
picture of above. A potential solution to this problem is to replace the string that re- 
sides in four dimensional flat space-time with a holographic string in non-flat higher 
dimensional space-time. An argument in favor of this scenario is the fact that the 
behavior of the Wilson line in an AdS^ space admits a j which is similar to the 
behavior of the quark anti-quark potential in short distances, namely E ~ where 
a{t) is the running coupling. This naturally calls for a holographic description in 
terms of a background that admits a conflning behavior [5] for long strings and which 
is asymptotically AdS. 

Precise lattice simulations of the quark anti-quark potential especially in the 
short distance range may serve as a test of the idea of weakly curved holography. 
Whereas for long strings the contributions of massive stringy modes are negligible, 
for short strings the impact of the massive modes becomes more important. In 
holographic models massive modes arise from fluctuations of the string along the 
radial direction, compact directions and along fermionic directions. Thus high pre- 
cision proflle for short distances will enable one to rule out or conflrm the idea of 
describing the Wilson line in terms of a string in weakly curved higher dimensional 
space-time. 
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The goal of this paper is to search for a holographic background that admits 
a Wilson line which is as close as possible to the potential deduced from lattice 
calculations and which at the same time is in accordance with other requirements. 
In [5] the conditions for a background to admit an area law behavior for long Wilson 
lines were stated. For short Wilson lines the holographic model is required to mimic 
the QCD perturbative behavior of the potential. To determine this condition is 
a complicated task. Here instead, in order to have a guide in the search of the 
appropriate background, we determine the conditions for a background for which 
the short distance potential behaves as j. 

In the UV the lattice simulates QCD in its perturbative regime. The holographic 
correspondence maps this regime to a highly curved gravitational background. Thus 
obviously the comparison we make to the holographic stringy result is on a solid 
ground only for large and maybe intermediate distances. One may wonder whether 
any reliable information can be extracted from the stringy picture in this case. Obvi- 
ously the background will be determined in this case by an action that includes higher 
order curvature terms and moreover the stringy calculation will include higher a' 
corrections. However, it is quite plausible that the maximally symmetric AdS back- 
ground will remain a solution thought the radius will differ from the classical value. 
For such a case, even though the N.G. is not an adequate framework of calculation, 
due to scale invariance the energy should scale as a |. 

It is well known since the seminal paper [6] that the classical string action that 
corresponds to a Wilson line diverges, and hence one has to invoke some renor- 
malization procedure. We discuss three different possible renormalization schemes: 
the holographic renormalization of [7], the Legender transform approach [8] and the 
mass subtraction [9]. It is shown that for backgrounds without an horizon^ the three 
schemes yield the same result but for cases that have horizon the last prescription 
differs from the two previous ones. From the large arsenal of holographic models we 
have analyzed the models of near extremal and branes [10], non-critical near 
extremal AcISq [11, 12] and the Klebanov Strassler [KS) background [13] ^. 

An outcome of this analysis is the fact that the quark anti-quark potential admits 
on top of the linear and j terms also a constant term. This constant term appears 
in all the holographic backgrounds we have analyzed and in all the renormalization 
schemes. It also follows from the lattice simulations, where it is defined as the dif- 
ference between the constants added to find the best fit with the simulation data in 
the large ^ range and in the short range. However, it seems that the constant term 
in the latter range is not well-defined in field theory as the theoretical perturbative 
expansion is not guaranteed to converge. Being aware of this point we suggest dif- 

^ Horizon in general relativity associates with the time component of the metric. Here and in 
the rest of the paper horizon will relate to a space coordinate rather than to the time. 

^Some similar aspects of the Wilson line in the KS background have already been discussed in 
[14]. 
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ferent ways to discuss the constant term in spite of this problem. The fact that the 
constants in the two different ranges are different is what renders this constant to be 
a physical quantity. A possible physical interpretation of this constant is in terms 
of masses attached to the string endpoints. A similar constant term has also been 
treated in the work of [15] on the zero point energy of renormalized Wilson loops, 
where the discussed action is a flat spacetime N.G. plus higher terms in covariant 
derivatives on the worldsheet. 

Quantitative comparison between the lattice simulations and the holographic re- 
sults shows that the Klebanov Stmssler background exhibits the smallest x square. 
However, the non-critical AdS^ model also admits a reasonable fit to the simulation 
results and it has the advantage of being asymptotically AdS, while the KS is only 
approximately asymptotically AdS and hence doesn't protected from quantum cor- 
rections in the UV . Both models are also invariant under four dimensional Poincare 
symmetry, which is an essential requirement. 

The paper is organized as follows. After this introduction, in section (2) we 
review the notion of Wilson loops in holographic models. We then discuss the con- 
ditions on the background to satisfy a coulombic potential at small distances. The 
issue of the holographic renormalization is discussed in section (4). We analyze sepa- 
rately the various renormalization schemes of [7] , [8] and [9] . A brief review of certain 
holographic models dual of Yang Mills like theories and discussion about their Wil- 
son line is presented in section (5). In section (6) prior to the performing of the 
comparison between the outcome of holography and the results of lattice simulation, 
we setup the units appropriate for the comparison. The latter is made in section (7). 
In section (8) we discuss the physical constant that follows from the simulations and 
the holographic models. We end the paper in section (9) with a brief summary and 
discussion of the findings of this paper and certain open questions. In particular we 
discuss the issue of the impact of quantum fluctuations on the holographic Wilson 
line. 

2 The holographic Wilson line 

We will investigate several holographic models, all of them are of the form of some 
high dimensional space with a boundary that consists some of the coordinates (from 
which only four should eventually be infinite), a radial coordiante that takes us from 
the bulk of the space to the boundary and can be viewed as energy scale and possibly 
some other coordinates transverse to the boundary's coordinates. Following [5] we 
will assume that the metric depends only on the radial coordinate such that its 
general form is 

ds^ = -Gtt {u) dt^ + Guu (u) du^ + ^^.11^.^1 {u) dx\ + G^^^^ {u) dx% (2.1) 
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where t is the time direction, u is the radial coordinate, X|| are the coordinates on the 
boundary and xt are the transverse coordiantes. We adopt the notation in which 
the radial coordinate is positive defined and the boundary is located at m = oo. In 
addition, an "horizon" may exist at u = Ma, such that spacetime is defined in the 
region ua < u < oo instead of < m < oo in the case where no horizon is present. In 
order to make our discussion most general, we will consider the case with an existing 
horizon. Later on, when we imply the general analysis to specific models, we would 
just set ua to zero in cases with no horizon. 

The construction we will examine is of an open string living in the bulk of the 
space with its both ends tied to the boundary. From the viewpoint of the field theory 
living on the boundary the endpoints of the string are the qq pair, so the energy of 
the string is related to the energy of the pair [5]. In order to calculate the energy of 
this configuration in the classical level we will use the notations and results of [5]. 
First, we define 

f{u) = Gu{u)G^^^^^^{u) ^^^^ 

g"^ (u) = Gtt (u) Guu (u) 
Upon choosing the worldsheet coordinates a = x {x is a. coordinate on the boundary 
pointing in the direction from one endpoint of the string to the other one) and r = t 
and assuming translation invariance along t, the Nambu-Goto action describing the 
string takes the form 

5"^^ = J dadr^ det [daXt'd,3X''G^^] 

dxdt^Jp (u (x)) + g^ (u (x)) (d^uf (2.3) 
T / dxC 



Then the conjugate momentum and the Hamiltonian are 

6C q'^{u)u' , 

V=^,= ; 2.4 

(n) + ^2 (n) n'2 

U = p-u' -C = -^-^ (2.5) 

As the Hamiltonian does not depend explicitly on x, its value is a constant of motion. 
We shall deal with the case in which u [x) is an even function, and therefore there is 
a minimal value uq = u (0) for which u' (0) = 0. At that point we see from (2.4) that 
p = 0. The constant of motion is therefore 

H = -f{u^) (2.6) 

from which we can extract the differential equation of the geodesic line 



du _^ f{u) {u)-P (no) 
dx g{u) f{uo) 
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and re-express the on-shell lagrangian (i.e. the lagrangian on the equation of motion) 
as a function of / (u) only 

Then the distance between the string's endpoints (or the distance between the 
"quarks") would be 



/{") 



where Uq is the minimal value in the radial direction to which the string reaches and 
Ub is the value of u on the boundary. The bare energy of the string is given by 

fdu\~^ ^ n , g(u) 
' ^ ' du- ^ ' 



Ebare (^o) = J dxC = J du ^ ^ j 



"0 / f(uo) Y 



V^(")/ (2.10) 



Ub 



f{uo)-i{uo) + 2 / dug{u)Jl- 



uo 



P K) 

P in) 



Generically, the bare energy diverges and hence a renormalization procedure is needed. 
Possible renormalization schemes are discussed in section (4). Here we follow [5] and 
use the mass subtraction scheme in which the bare masses of the quarks are sub- 
tracted from the bare energy. The bare quark mass is viewed as a straight string with 
a constant value of x, stretching from u = (or u = ua if there exists an horizon at 
Ua) to u = Ub, such that it is given by 

PUb 

rriq = / dug{u) (2-11) 

Then the renormalized energy would be given by 

^(wo) = /(wo)-n«o)-2/C(Mo) (2.12) 

where /C {uq) is 

/C(no)=/ dug{u)\l-.ll-'—^\+l dug{u) (2.13) 




In section (4) we make a discussion about the differences between different renormal- 
ization schemes. We emphasize that this result for the energy is only at the classical 
level and do not include quantum corrections. 

In order to reproduce the QCD heavy quarks potential we first have to demand 
the holographic models to reproduce the asymptotic forms of the potential. This 
leads to several restrictions on the forms of the / and g functions. The condition for 
confining behavior at large distances was derived in [5]: 
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Figure 1. The physical picture arising from the present construction. Left: The confining 
hmit is approached as the string is close to the horizon, then most of the string lies on 
the horizon and implies a string-like interaction between the quarks (the endpoints of the 
string). Right: The conformal limit is reconstructed when the string is far away from the 
horizon such that it is ruled by the geometry near the boundary. 

1. / has a minimum at Umin and / {Umin) 7^ or 

2. g diverges at Udiv and / {udiv) 7^ 

Then the string tension is given by / {umin) or / (udiv), correspondingly. The second 
asymptotic is perturbative QCD at small distances. In section (3) we derive the 
conditions on the background to reproduce the leading perturbative behavior of 
QCD, which is Coulomb-like (throughout this paper we refer to a | behavior as 
Coulombic). 

The physical picture arising from this construction is as follows. The confining 
limit is approached as uq — )■ ua, then most of the string lies on the horizon^ and 
implies a string-like interaction between the "quarks" (the endpoints of the string) 
with a string tension 

a = fM (2.14) 

and a linear potential 

E = a-i-2K (2.15) 

where k = IC (ua) is a finite constant. On the other hand, the conformal limit is 
approached as Uq/ua — 00, then the whole string is far away from the horizon and it 
is ruled by the geometry near the boundary. In other words, this limit correspondes 
in some sense to taking u\ to zero, such that the scale parameter disappears and the 
conformal theory is reconstructed. 

^For simplicity we depict the picture only for the case where horizon is present, but a similar 
picture arise also in cases with no horizon. For more extensive discussion see [16]. 



- 6 - 



3 The conditions on the background to satisfy a coulombic 
potential at small distances 



In order to generate the QCD heavy quarks potential in an holographic model it 
should first reconstruct the potential at asymptotics. The conditions on the holo- 
graphic background to reproduce the confining asymptotic at large distances were 
derived in [5] . As was discussed in the introduction, the holographic dual of the per- 
turbative region is necessarily highly curved and hence unreliable. Nevertheless we 
would like to derive the conditions on the background to reproduce at small distances 
a potential of the form 

(3.1) 

The full derivation is detailed in appendix (A). 

The first step would be to expand the / and g functions around the boundary, 
which we will take to be at infinity Uf, = oo, 

/(n) = afcn^- + 0(n^-i) 

g{u) = g (oo) + hju-^ + O (3.2) 
fc,j>0 

We assume here that / {u) diverges on the boundary as a power law and that g (u) 
tends to a finite value on the boundary. Then we can calculate the length of the 
string as a function of its minimal radial position, uq, using (2.9). Since we are 
interested in the small distance behavior, i.e. Uq which is close to the boundary, the 
calculation will be done in the region Uq ^ Ma, there we find 

i (no) = (^C^2M. (oo)) + O (-^) (3.3) 

V flfc / (no) VK) / 



where 



1 



C«.M^/^ (3.4, 



is a convergent integral. Next we calculate K, (uq) in the same manner 



K. (wq) = uog (oo) 



1 + lc!,,^2k M 



+ O W"^) (3.5) 



(where C'g;. 2fc (oo) is another convergent integral defined in appendix (A)) and use 
(3.3) to find the inverse relation of Uq as a function of £ 



1 

-1 



Mo = C2k,2k (oo) - (3.6) 
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in order to express AT as a function of the length 

1 



2fc,2/t 



OO 



[2C2k,2k{oo)]—[g{oo)]— { — 



(3.7) 



Then we are ready to substitute these results into the expression for the energy of 
the string (2.12) as a function of its length to get 

Eii) = fiuo (£))■£- 2/C {£) 



1 



(3.8) 



where is a /c-dependent dimensionless number 



iV, 



1 + 2C2k,2k (oo) 



2fc,2fc 



OO 



[2C2fc,2fc (oo)]' 



(3.9) 



The bottom line is that the dependence of the energy on the length in the region 
close to the boundary, i.e. small distances, is 



E 



1 \ ''^^ 



(3.10) 



We conclude that in order for the potential to satisfy a coulombic form near the 
boundary, one has to demand k = 2. In other words, the condition on the backround 
is 



/(n>nA)~n2 (3.11) 

We emphasize that this analysis is valid as long as / diverges on the boundary as a 
power law and g is finite there. An example for a case that fulfills this ansatz is the 
famous AdS^ x 5*5 solution in which 

^^"^ = 2^(^ 
^^"^ = 2^ 

4 Renormalization 

Generically, the bare energy of the string in the holographic construction (2.10) is 
divergent. From the bulk perspective it is a result of the infiniteness of the boundary. 
The string's endpoints are tied to the boundary and the energetically preferred state 
is when the string extends to the bulk^. Since the boundary is radialy located at 

^There may be gravitational solutions in which this is not the case and the energetically preferred 
state is when the whole string extends to the boundary only, but then the Wilson loop is just the 
flat space Wilson loop. We do not consider such a case as we are interested in holographic Wilson 
loops. 



1 9^1 
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infinity, this would mean that the string's length is infinite, and this is the origin of 
the divergence. From the field theory living on the boundary perspective the infinite 
bare mass of the quarks is the source for the divergence. The most simple and naive 
approach would be to subtract this bare mass [6], as one can think of the quark as 
a straight string stretched from u = (or m = ma if there exist an horizon at u\) to 
u = such that its mass is given by 



du g [u) 



(4.1) 



MA 



That approach was taken in (2) following [5]. However, there may be other regular- 
ization schemes as well, and we want to discuss and compare them. In this section we 
examine two other schemes: the holographic scheme [7] and the Legendre transform 
scheme [8]. We calculate analytically the subtracted part of the energy in each of 
the schemes, and then discuss the results. 

The calculation is done using the expression for the bare energy (2.10) and 
assuming the following assumptions on the / and g functions: 

1. / is a monotonic increasing function 

2. is a constant function or a monotonic decreasing function which does not 
diverge on the boundary 

The first one is the same as the ansatz for the / function in [5], while the second 
ansatz is more restrictive, but in fact is not needed in our analysis of the Legendre 
transform scheme. 

4.1 Holographic renormaliztion 

The idea behind this scheme [7] is to expand the bare energy around the boundary, 
identify the infinite terms and subtract them from the original expression of the 
bare energy. Essentialy we remove infinities resulting from the infiniteness of the 
boundary. Therefore we will not lose any information about the subtracted part if 
we will consider the bare energy with uq which is far from the horizon. Then g (u) 
is approximately constant in all the range of integration 



rut ^ 
Ebare (%) = "^9 (Ub) / du 



(4.2) 



Now let us divide the integration into two parts 



Ebare [Uq) = '^9 (Ub) 



du 



+ 



/(«) 



f(no) 

/(«) 



du 



Mo 



1 



f{n) 



(4.3) 



^In this section we will keep the notation in which the boundary is located at w = Ub- Later on 
we will set Uh = oo. 
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such that / (n) ^ / (uq) above so we can simphfy this expression 



Ebare {Uq) = 2g (Ufe) 



du 



du 



1 



no 



1 - 



2ukg (ub) - 2g (ub) Uq 



u 



/ dy 

Uq Ji 



1 



f /(»()) 
V /(«) 



(4.4) 



where in the last transition we have changed the integration variable to y = The 
second term does not contains boundary infinities, because its only dependence on 
the boundary is via g (ub) which is finite in our ansatz. So even if the second term 
includes divergent parts, this scheme will not take care of it. Writing the bare energy 
in that form we isolate the infinity resulting from the boundary such that according 
to the prescription of this scheme, the subtracted part of the energy would be 

Esub = 2ubg (ub) (4.5) 
4.2 The Legendre transform scheme 

The idea behind this scheme [8] is to consider the legendre transform of the energy 
as the renormalized energy 



E, 



renormalized 



J dx {C — p ■ u) 



such that the subtracted part of the energy is 

Esub = J dx{p- u') 
Now, since p does not depend explicitly on x (see (2.4)) we can write it as 



sub 



(4.6) 



(4.7) 



(4.8) 



which means that the subtracted part of the energy is a boundary term, and therefore 
we can expect to have a similar result as in the holographic scheme. Changing 
variables from x to u {dx-^ = du-^) and substitute the momentum (2.4) on the 
geodesic line (2.7), we have 



Esub = du— 

Juo 



du 



2ug (u) 



2ubg (ub) \ 1 



(4.9) 



2ubg (ub) 
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which indeed result as the holographic scheme. Pay special attention to the fact that 
we have not used here the ansatz for the g function stated in the beginning of this 
section and still got the same result as for the holographic scheme. 

4.3 The mass subtraction scheme 

This scheme suggests to simply subtract the bare mass of two free quarks, each of 
them is represented by a straight string stretched from the boundary to the horizon. 
Therefore the subtracted energy is given by [5] 



E 



sub 



du g [u) 



UA 



Let us seperate the integration into two regions 



sub 



du g (m) 



du g (m) 



(4.10) 



(4.11) 



UA 



where above u^, g is approximately stabilized on its boundary value 

"Ut 



Esub = '2g (uh) / du 

J Um 

= 2ubg [ubj - Ua 



du g (u) 



UA 



Ur 



2—g{ub)-2 



(4.12) 



dygiy) 

where in the last transition we have changed the integration variable to y = The 
first term is essentialy the same one as of the two other schemes, and is a result of 
the infiniteness of the boundary. The second term is a feature of the horizon and it is 
finite since it does not contain boundary infinities (because we isolated them) and its 
only dependence is on the horizon characteristics, which cannot cause a divergence of 
the energy (because the curvature is not singular). Its exact value can be calculated 
for a given model. We conclude that the subtracted energy in this scheme takes the 
form 

Esub = 2ubg (ub) - A ■ Ua (4-13) 
where A is a constant finite number that depends on the exact form of the horizon 



A = limi 



2^giub)-2 
Ua 



Um 

UA 



dygiy) 



(4.14) 



In the non critical AdSn+i solutions that will be discussed later, for example, one 
finds 



A{n) 



(4.15) 
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U = Ub 



u = 




£ Ub — UA ■ function(horizon) 




Figure 2. Left: In the case in which the renormahzation scheme ignores the horizon, one 
can think of this scheme as of subtracting the length of a string streched from the boundary 
to the end of space at u = 0, which is proportional to the location of the horizon, Uf,. 
Right: In the mass subtraction scheme one subtracts the length of a string streched from 
the boundary to the horizon, which is shorter by a finite amount from the string of the 
first the effect of the horizon must be finite. On dimensional grounds, this finite 

difference must be proportional to ua and it depends only on the characteristics of the 
horizon [function (horizon) is a finite dimensionless number). 

4.4 Conclusions 

We see that the mass subtraction scheme differs from the two other schemes by 
a constant term in the expression for the subtracted energy. The fact that differ- 
ent schemes differ only by a finite term is true in field theory and here we have a 
geometrical realization for that feature, as explained below. 

The holographic and Legendre schemes are taking care only of the infintiness of 
the boundary, and do not take into account the existence of the horizon, while the 
mass subtraction scheme does take it into account. Indeed, in cases with no horizon 
we set ua to zero (according to the prescription in section (2)) and the difference 
between the schemes disappears. As was mentioned in the beginning of this section, 
the divergence is a result of the infiniteness of the boundary. The difference between 
the schemes is only the horizon, which cannot cause any divergences of the energy 
since the curvature is not singular. Therefore we conclude that the difference between 
the schemes must be realized as a finite affect. 

An important implication of this conclusion is that the analysis of the behavior 
of the Wilson loop at asymptotic regions (in [5] for the confining region and in section 
(3) for the coulombic region) is independent of the regularization scheme (as far as 
those three schemes are considered). 

Furthermore, out of the fact that the UV region is eventually decoupled from the 
IR region as we take the limit Uf, — )■ oo, one can restrict the form of the subtracted 
parts. Each of those regions has its own scale parameter and u\, respectively) 
and therefore properties of each of them would depend only on the relevant scale. 
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The divergent part of the energy is a feature of the UV region (since the divergency 
is a result only of the infiniteness of the boundary) such that on dimensional grounds 
it would be proportional to the location of the boundary, Ub- On the other hand, the 
finite difference between the schemes would be proportional to u\ since it is a result 
of the IR properties of the background. We describe this picture in figure (2). 

5 The Wilson line of the holographic models 

So far we have made a general discussion about the Wilson line in holographic back- 
grounds. As we aim to compare the QCD Wilson line to the one of holographic 
theories, we have to consider specific models. In this section we briefly review few 
representative holographic models, which we latter compare to lattice simulations, 
and discuss their Wilson line. Among the large class of holographic models we have 
choose to consider the near extremal and D4 branes, non critical AcISq and the 
Klebanov Stmssler construction. The reason for this choice is that these backgrounds 
admit a Wilson line which is in principle close to the QCD Wilson line (at least on 
asymptotics) , as we will discuss below. In addition, they are all solutions of string 
theory or of supergravity. Other examples have been analyzed in the literature, 
among them also phenomenological models and backgrounds which are not solutions 
of a particular action. In this paper we choose not to deal with such cases since one 
can always engineer a background which will admit some required features without 
being a solution of an action. Instead, we want to ask what is the solution that 
admits a Wilson line which is closest to the QCD Wilson line. Nevertheless, for a 
matter of completeness, we shortly review several other results; The Wilson line of 
the phenomenological Holographic QCD [17] has been analyzed in [18] where a good 
agreement to the lattice results was found. In [19], the Wilson line of some hard- wall 
scenario was calculated and it was found to admit a Cornell potential form®. 

Near extremal Dp brane 

In the supergravity background associated to near extremal Dp brane [20] one finds 

^ ^ ' \27ra' J 1 _ [^Y'^ 

These backgrounds satisfy the condition for confinement for all values of p. On the 
other hand, using the analysis of (3), we see that the potential in the small distances 
region takes the form 

-Gf 

^We thank Andrew Zayakin for referring us to this paper. 
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Therefore, only for p = 3 the potential behaves as | in the small distances region. In 
that case there are four coordinates on the boundary, where one of them is compact- 
ified. For that reason, on energy scales smaller than one over the compactification 
radius 

E < -s (5.2) 

-'^compactification 

the theory on the boundary is effectively three dimensional. The conclusion is that 
the field theory on the boundary of the near extremal brane is four dimensional 
at small lengths scale (or equivalently high energies scale), where the potential is 
Coulombic, but it is effectively three dimensional on large lengths scale (low energies), 
where the potential is confining. In order to finish with an effectively four dimensional 
boundary at large distances, one needs to consider the near extremal D4 brane, 
following the same argument. In that case the boundary theory on small length 
scales would be five dimensional and the potential there will take the form ~ ^. 

Non critical AdS solutions 

The metric of a non-critical AdSn+i solution of type II string theory [11, 12] with 
one of the spatial coordinates compactified on a circle is given by 



-df + dxi + {1- { — ] ]dx 

u 



comp 



1 / RAds \ du 



2Txa' \ u y (1 - (v)") 



2Txa' \RAds, 
i = l,...,{n-2) 

(5.3) 

where u\ is the energy density on the brane. Therefore the / and g functions take 
the following forms 

I f u 



fin) 



^AdS 



^ (5.4) 



This background satisfies both conditions for confinement at long distances and j 
potential at short distances, for any value of n (larger than 2). For these conditions 
we have to add the demand that the effective boundary theory at long distances 
would be four dimensional. Therefore, a good candidate would be the non-critical 
AdS^ which has a five dimensional boundary with one compactified coordinate (and 
hence a four dimensional effective field theory on the boundary, following a similar 
argument as in the D branes discussion). 

Klebanov Strassler 

The supergravity solution of the deformed conifold is of the following form [13, 21] 

ds^ = (r) dxl-^23 + (''") (5.5) 
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where dsi is the metric of the deformed conifold 



dsi 



A/3 



-K(t) 



3K^ (t) 



cLt'^ + [g^Y] + cosh 



+ sinh 



e is the energy scale and the functions h (r) and K (r) are given by 

h{T) = {gMa'f^'/h-^/^I (r) 
1 



It) = 



°° X coth X 
dx- 



Kir) 



sinh X 
(sinh {2t] 



(sinh 2x — 2x) ^ 



T 

2 

(5.6) 

(5.7) 
(5.8) 

(5.9) 
2 3 sinh r 

r is the radial coordinate, it is a pure numer (i.e. unitless) running from zero to 
infinity on the boundary. Near the boundary, i.e. at large r, the h (r) and K (r) 
function takes the form 

/i (r > 1) = {g.Ma'f 3 ■ 2^/^e~^/^Te-'s^ (5.10) 

K (t :$> 1) = 2h-^^ (5.11) 
From the form of the metric we conclude that the / and g functions are given by 



2\2 



(^^) +(/) 



1 



(5.12) 



Following [21] one can check that the KS background is confining since / (r) 
has a minimum at r = and it is finite there / (0) 7^ 0. On the other hand, near 
the boundary (r ^ 0) this backgroud do not follows the form (3.2) required for the 
analysis in section (3). However, it is approximately AdS in this region, therefore we 
expect it to generate a potential which is approximately j at short distances. Near 
the boundary the KS metric takes the form 



T^dr 



(5.13) 



The behavior of the wrap factors is dominated by the exponential term. Therefore, 
by neglecting the power law terms of the wrap factors we change the geometry just 
by little. The resultant approximate geometry is then 



ds' ~ e'/^'el^dxl,,, + dr^ 



(5.14) 



which can be brought, by an appropriate coordinate redefinition, to the form of an 
AdS geometry 



ds' 



u \ 



RJ 



du^ + 



(5.15) 
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Hence, since the AdS geometry generates a | potential, we expect the KS background 
to generate a potential which is approximately j near the boundary. Then the 
arising question is how do the corrections to the approximate j potential look like? 
To answer this question we use the theorem reviewed in (2) (the exact analysis is 
detailed in appendix (G)). Plugging the discussed geometry (5.13) into the expression 
for the length of the string (2.9) one finds 



,2/3. 



9sMa' 1 
2/3 3 



/ (To) 

In the same way the Kappa function (2.13) takes the form 

JC (To) ~ e2/3e^« 
Then the energy of the string (2.12) would be 

E (ro) ~ -e^/^el-o 



(5.16) 



(5.17) 



(5.1^ 



In order to derive the dependence of the energy on the length we need to invert 
the function i (tq) into tq {£) and plug it into the expression for the energy E {£) ~ 
_g2/3g3To(£)_ ^]-^ig jg very difficult to perform. Instead we will expand i (tq) around 
large tq and then invert the relation. For that it would be usefull to take the logarithm 
of the length 



In 



The leading order 



In 



would result in the leading ^ behavior 



E 



(5.19) 
(5.20) 

(5.21) 



as expected. In order to calculate the leading correction we would plug the leading 
order relation (5.20) into the subleading term of the exact expression (5.19) 



In 



{9sMa') 



--To + In y^- In [{gMa'V e^] 



(5.22) 



to get the dependence of Tq on 
1 



3 " 



In 



{gsMa'Y^ e^'H + In ^- In [(^.Ma')"' e^'H] (5.23) 
Inserting this into the expression for the energy we finish with a potential of the form 



E{i) 



In {gsMa') 



-1 e2/3£ 



(5.24) 



In principle, there could appear also a constant term, as will be discussed in section 
(8). In appendix (G) we explicitly calculate it. 
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6 Units of measurement 

In order to make a quantitative comparison between different models one has to 
discuss the issue of units of measurement. In this section we discuss two issues 
concerning this subject. We start with an abstract discussion and then specialize to 
the models under investigation. 

The first issue is the question of how to determine the system of units when 
infinitely distant scales are introduced. In a theory with only one scale parameter 
(say, the fundamental string length 1^) every physical quantity Q will be expressed 
as some power of the only scale 

Q ~ {IsT (6.1) 

Now consider the case in which another scale R is introduced. When the two scales 
are very far apart Q would have to be expressed as some powers of both 

Q ~ {RT (6.2) 

If we further seperate the two scales by an infinite gap, there would be a unique 
choice of the powers pi and p2 such as to get finite answers, and this is the way one 
has to choose his system of units. Any other system would make no sense since it 
would result with infinite values for physical quantities. 

The second question is how to compare different models which have different 
scale parameters. Each model has its own scale parameters and therefore the im- 
mediate system of units is different in each of them. In order to quantitatively 
compare different models one has to describe the quantities under comparison in a 
model-independent way, i.e. they have to be expressed in terms of universal physical 
quantities that are present in all the models under investigation. 

Specializing to the framework we explore in this paper, the physical quantities we 
are interested in are the energy of the string and the distance between its endpoints 
on the boundary, the scales are the fundamental string length Ig and the radius of 
the space Rspace, and the fact that the two scales are infinitely far away from each 
other is exactly the statement of weakly curved holography 

a' = I', « (6.3) 

It is clear now that since Ig and Rgpace do not exist in the lexicon of the gauge theory, 
we need to express the energy and length in terms of universal physical quantities 
which are present in both frameworks, stringy holography and gauge theory. Since 
we are dealing with aspects of strong interaction, the natural universal physical 
quantities are the string tension and the glueball mass (because they must exist in 
any model that aims to describe strong interactions). In this language, the SUGRA 
limit is 

^ > ^llueball (6.4) 
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Being expressed in such a way, we would be able to compare the Wilson line of these 
two frameworks on equal footing. 

As a realization of this general discussion, we first consider the near extremal 
Dp brane. In appendix (C) we apply the theorem presented in section (2) to find 
the expressions for the energy and length of the string in this case. We find that the 
energy of a classical string is well-defined when it is measured in units of 

[^1 - ^ (6.5) 

where ma is the scale parameter and a' is related to the fundamental string length. 
The length of a classical string is well-defined when it is measured in units of 

M = u, (6.6) 

where RacLS is the radius of the AdS space and is of the same order as ujy. In this 
system of units, following equation (2.14), the string tension will take the value 

1 

a = — 

2ti 

in units of (the units of the string tension are essentially the units of the energy over 
the units of length) 

[E] If UA y-^ 

[i] a' \RAds) 

Any other choice of system would result in an infinite value for the string tension in 
the SUGRA limit. The glueball mass scale is related to the radius of the compactified 
coordinate ^ 

T^gb ~ (6.7) 

■^compactified 

which is uniquely determined from the demand of avoiding a conical singularity. 
In appendix (B) we review this calculation and show that in the case of the near 
extremal Dp brane background it is given by 

7-p 

47r / Ma \ ^ 



27rR^ ^ = UA (6.8) 

compacti J tea T ^ \ ID I ^'^ 

' — p \^AdS/ 

Then the glueball mass is 

Ua \riAdS / 

7-p 

The relation between the string tension (o" = {^^EZw) ^ glueball mass 

would take the form 

a' 

rrigb ~ — a (6.10) 

Ua 
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Now we can express the energy and length in terms of a and mgf, 



[E] = ^-^ (6.11) 
a nigb 



M = ( 1^ ) ~ — (6.12) 



ua y mgb 

Being expressed in such a way, the energy and length are well-defined and their values 
on both sides (holography and gauge theory) are compareable (given that they are 
expressed in a similar way in the gauge theory). The exact translation is derived in 
the next section. 

For the case of non-critical AdSn+i a similar analysis holds. The string tension 
is given by 

-^(^V (6.13) 



27ra' V R 



l-AdS 



.2 



AdS 



n UA 



As the periodicity of the compactified dimension is given by /3 = 27rR 
(see appendix (B)), the glueball mass would be 

^9b ~ (6.14) 

^AdS 

and the same relation (6.10) between rrigb and a holds. Using the appropriate units 
for the energy and length found in appendix (E) we reconstruct (6.11) and (6.12) for 
the present case 

[E] = ^-— (6.15) 
a rrigt 

[i^ = ^r^J_ (6.16) 
Ma mgb 

Next we consider the Klebanov Strassler construction. The string tension is given 
by (2.14) 

4/3 4/3 

a = / (0) = /."i (0) = 2-1/3/ (0) - 0.57^-— (6.17) 

and the glueball mass is [13] 

£2/3 

where e is the scale parameter and Qs is the fundamental string coupling and M is 
the rank of the symmetry group. Using the form of the metric (5.5) one can find the 
appropriate units for the energy and length of the string^ and express them in terms 
of the string tension and the glueball mass 



^ Since the radial coordinate is dimensionless, the units of the energy are essentiahy the units of 
the g function whereas the units of the length are the units of the g function over the units of the 
/ function. 
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a 


TD I, = T77r»-4--^ 


^3 brane 


1 ( UA ^ 
27ra' \^ RAdS J 


'2 


Vll.59-J^ = 3.4-^ 

-^AdS ^AdS 


D4 brane 


1 ( UA \ 

2-Ka' \RAdS ) 


2 


,,1/2 ,,1/2 

(9-39 ■ 1^) ^ = 2.24^ 


AdSe 


I ( UA ^ 
2-Ka' \RAdS ) 


2 


(6.34-1^) ^2.52^ 

^ '^''^ ^AdS ^AdS 


KS 


0.57 


0.86 



Table 1. String tension and lowest 0''""*" glueball mass of the various holographic models. 



[E] = e^/^r.— (6.19) 
rrigb 

gsMa' 1 

The exact translations between the natural system of units of each holographic 
model to the universal one is derived in the next section. 



7 Comparison with lattice gauge theories 

In this section we will compare the potential of the holographic models to lattice 
results of SU (2) and SU (3) gauge theories without fermions [22, 23]. In order to 
do so we first have to express everything in terms of universal physical quantities, as 
described in section (6). Our next task would be to determine the exact translation 
between the natural system of units of each holographic model to the universal 
one. In table (1) we summarize the results for the string tension and the lowest 
0"'"'" glueball mass^ of the various holographic models. The lowest 0"'""'' glueball 
mass for the D3, D4, AdS^ and KS models was calculated in [24], [25], [12] and 
[26], respectively. The units of [26] 's results are explicitly resolved in [27]. Then it 
is simple to invert the relations such as to result with a dictionary of translation 
between the natural system of units of each holographic model to the universal one, 
see table (2). 



The next step would be to express the system of units used by the lattice con- 
ventions in terms of the universal physical quantities discussed above. In the lattice 

*The reason for choosing the 0^^ state for the universal system of units is that in SU (N) gauge 
theory it has the lowest mass. 
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[E] 




brane 


^ = (3.4 ■ 27r) ^ 


1 

^AdS _ 2 1 1 


D4 brane 




f ,f = 2.24 1 


AdS^ 


^ = (2.52 • 27r) ^ 


^ids ^ 2.52 ^ 


KS 


e2/3 = 1.51^ 


= 0.86 1 



Table 2. Dictionary of translation between the natural system of units of each holographic 
model to the universal one. 



conventions, length is measured in units ro related to the string tension via the next 
expression [23] 

arl = 1.65 - ^ (7.1) 
The value of the lowest 0'^"'" glueball mass in D = 3 + 1 dimensions is given by [28] 

mgiuebaii = mo++ = (In ■ y/a (7.2) 



d 



N 



3.87 N = 2 
3.65 iV = 3 



Hence we can deduce the relation between the original lattice length's unit to the 
desired ones 

ro = dnJl-Qh-^— (7.3) 

In the same manner the translation of the energy units, which are given in the lattice 
convention by ^, is done using the next expression 

1 a 



ro ^1.65 - nigb 

The exact expressions for the potentials in Dp brane, non-critical AdSn+i and 
Klebanov-Strassler backgrounds are derived explicilty in appendices (C), (E) and 
(G), respectively. They are presented together with the lattice data in terms of the 
universal system of units and using two different conventions 

1. Fixing the zero energy point (as in the conventions of [23]). 

2. Comparing the confining asymptotic behavior of all the graphs, 
in figures (3) and (4), respectively. 
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Heavy Quarks Potential - Fixed Zero Energy Point 




Figure 3. The heavy quarks potential of the holographic models and the lattice data 
expressed in the appropriate physical units and presented using the fixed zero energy point 
convention. 
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Heavy Quarks Potential - Unified Confining Asymptotic Behavior 




Figure 4. The heavy quarks potential of the holographic models and the lattice data 
expressed in the appropriate physical units and presented using the convention in which we 
compare the confining asymptotic behavior of all the graphs. 
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Convention— > 


Fixed zero 
energy point 


Best fit 


Unified confining 
asymptotic behavior 


near extremal brane 


24.6 


24.4 


34.8 


near extremal D4, brane 


108 


108 


120 


non critical AdS^ 


4.14 


3.98 


5.89 


Klebanov Strassler 


1.4 


1.33 


14.5 



Table 3. values for the measure of correspondence between the hofographic models' 
potential and the lattice measurements for SU(3). The calculation was done using three 
different conventions, as detailed in the text. 

In order to quantify the measure of correspondence of the holographic models' 
potential to the lattice data we will use Pearson's test 

X' = ^lj^^^^ (7.5) 

where n is the number of the lattice data points. The computation is done using 
the two conventions described above and in addition we present another convention 
in which the height of the theoretical potential curve is adjusted to result with a 
minimal i-S- the best fit possible. We calculate only for the correspondence 
with the SU(3) simulation^. The results are presented in table (3). 



From the results we can determine that two of the models we have examined 
correspond significantly better to the lattice measurement of the Wilson line than 
the two other ones. The good models are the Klebanov Strassler and the non- 
critical AdS^ backgrounds which have significantly much better values than the 
near extremal D branes. On the other hand, the D4 brane background generates a 
potential which do not fit at all to the lattice potential, with x^'s greater than 100 
for all the conventions of comparison we have used. 

Out of this summary we can already conclude that the requirement for a Coulom- 
bic behavior at short distances proves itself in the sense that the models which admit 
at short distances a potential which is either exactly or approximately of the form 
~ I are also in a good agreement with the lattice measurements, compared to the 

brane which admits a ~ potential in this regime and has a value which is 
far above all the other. Hence it makes sense to use this requirement when searching 
for an holographic model which mimics QCD's Wilson line. 

The best fit convention of comparison gives almost identical results to the fixed 
zero energy point convention. In both, the model which has the best correspondence 

^We read [22] 's results for the SU(2) simulation out of the graph since they are not given 
explicitly, therefore they will be used only for a qualitative picture and not for quantitive measures 
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with the lattice results is the KS with a remarkable x square value, but nevertheless 
the AdS^ also fits quite well. In the third convention {unified confining asymptotic 
behavior) they flip over as the ^^5* becomes secondary to the AdS^ which best fits. 

8 The constant terms 

At the beginning of this paper we have made a statement that at large distances the 
potential is linear and at small distances it is Coulombic. However, there is still a 
possibility of an additional constant term for each of these asymptotic behaviors, i.e. 

Eiong = Cjl + CL (8.1) 

Eshort ^ -j + Cs (8.2) 

In a sense, these constant terms characterize the behavior of the potential at interme- 
diate distances. Therefore, a comparison between such constants of the holographic 
models to the constants of the gauge theory would constitute a measure for the 
similarity of the potentials at intermediates. 

In principle, each of the constants [Cl and Cs) for its own is meaningless and 
only their difference 

AC = Cl- Cs (8.3) 

is meaningful, since the energy can always be shifted by a constant. Therefore one 
should expect to consider only their difference. However, since we have fixed the zero 
energy point there is still a possibility to look at each of the constants separately. 
Of course, that would depend on the way we have fixed the energy, but we will still 
find it usefuU to compare the large distance constant of the different models and the 
lattice, as will be described later. 

Next, we will check whether such constants do appear in the holographic models. 
In appendices (D),(F) and (G) we calculate the values of and Cs for the different 
cases, before fixing the zero energy point as in [23] 's notations. The results we get 
after fixing the zero energy point are listed in table (4). Of course, the results for 
AC do not depend on how we have fixed the zero point energy, therefore they are 
explained by the mentioned appendices. 



From these results we can already conclude about the constant term in hologra- 
phy by itself. The values of AC in all the models under investigation are all non-zero, 
negative and very close to each other (AC ^ —5^^), such that several remarks are 
to be made: 
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Cl " 


Cs " 


AC " 

mgt 


brane 


0.59 


7.39 


-6.8 


brane 


1.35 


5.3 


-4.48 


non critical AdS^ 


-1.32 


4.02 


-5.34 


Klebanov Strassler 


-1.44 


2.93 


-4.37 



Table 4. The values of the constant terms for the different models in the fixed zero energy 
point convention. 

• It seems that generically holography generates such a constant term. We can 
get some intuition about this feature of holography by thinking of how one 
would generate it in a flat space-time Nambu-Goto theory. While at short 
distances a tachyonic instability takes over such that there is no sense in talking 
about a constant term in this region, at long distances it does make sense, but 
as happens - it does not shows up (see eq. (1.1)). The only way of generating 
such a term at long distances is by adding a boundary term to the Nambu-Goto 
action. Therefore, since holoraphic models have a boundary by definition, we 
can expect them to have a constant term of the discussed type. 

• We interpret (—AC) as the physical mass of the quarks which remains after 
subtracting their bare mass. Therefore it is crucial that the values of AC in 
all the models have the same sign. 

• The fact that all the holographic models admits a similar AC value is very 
interesting. 

Now, we will derive the values of the constant terms for the gauge theory, using the 
lattice results. By adjusting the theoretical prediction 

E = ai-^-^ + CL (8.4) 

(where the first term is just the linear potential and the second one is the well-known 
Luscher term) to the measurements, we find Cl- Finding Cs appears to be a subtle 
issue. In principle, one should compare the theoretical prediction of perturbation 
theory to the measurements, and derive Cs in the same manner as for Cl. However, 
as was found in [29], the pertubative expression of the potential (up to three loops) 
fails in describing the potential even at the shortest lengths reached by the lattice 
simulations. Therefore, by adjusting this theoretical prediction to the shortest lattice 
measurement, we can only get a rough estimate for C5. We are required to this point 
in the next paragraph, but first let us present the results, listed in table (5). We use 
theoretical predictions of 0,1 and 2 loop orders of perturbation theory, in order to 
see how much Cl changes as we go to high order expressions. 
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Cl " 


Cs " 


AC " 


Lattice SU(3) 0-loop 


-0.63 


-2.24 


1.61 


Lattice SU(3) 1-loop 


-0.63 


6.31 


-6.94 


Lattice SU(3) 2-loop 


-0.63 


2.48 


-3.11 


Lattice SU(2) 0-loop 


-0.43 


-0.01 


-0.42 


Lattice SU(2) 1-loop 


-0.43 


11.67 


-12.1 


Lattice SU(2) 2-loop 


-0.43 


5.94 


-6.37 



Table 5. The values of the constant terms for SU(2) and SU(3) gauge theories. We used 
the data of [22, 23] to calculate them. 




3 Number of Loops 



Figure 5. AC values as were calculated from comparison of the SU(3) lattice measurements 
to the theoretical predictions of 0,1 and 2 loops orders of perturbation theory. 



As was discussed above, the process of finding the constant term at short dis- 
tances (and therefore also AC) out of the lattice measurements is a subtle issue. One 
way of resolving this issue is to consider only the constant at long distances, which is 
well defined since we have fixed the zero energy point, but also somewhat arbitrary 
defined, since this fixing is arbitrary. By examining the values of the constant terms 
at long distances one can notice that all the holographic models has Cl's which are 
of the same order as the one extracted out of the lattice data for the gauge theo- 
ries. Moreover, our two favored candidates {KS & AcISq) have negative Cl as in 
the gauge theories (while the D branes have positive Cl), which may be regarded as 
another support in them. We emphasize that this conclusion should be taken with 
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a grain of salt since it depends on how we have fixed the zero energy point, but it is 
still interesting. The other approach we take here is to look on the values of AC as 
were calculated from adjusting the theoretical predictions of 0,1 and 2 loops orders 
of perturbation theory to the lattice measurements and try to estimate to which 
value it would converge if we were able to use the theoretical expression of higher 
orders. From the results of the SU(3) measurements (plotted in figure (5)) the guess 
is that AC would converge to a value somewhat around AC = —5-^. Considering 

'^gb 

this value as a basis for comparison, all the holographic models under investigation 
definitely have a very similar AC value. 

9 Concluding remarks & future directions 

Following the statement made in the introduction, the goal of this paper is to search 
for a holographic background that admits a Wilson line which is as close as possible 
to the potential deduced from lattice calculations and which at the same time is 
in accordance with other requirements: Confinement at long distances, Coulombic 
potential at short distances (by which we mean j behavior) and invariance under 
four dimensional Poincare symmetry. 

We summarize by suggesting the Klebanov Strassler background as the holo- 
graphic model which best fits the lattice measurements and in the same time is in 
accordance with the other requirements we have mentioned above^^. Nevertheless, 
the non-critical AdS^ also has an appreciable correspondence and it also admits the 
other requirements, so we cannot rule it out. However, from this point of view, the 
near extremal D branes can undoubtedly be excluded. 

We turn now to final remarks and future directions. The first remark concerns 
the rank of the symmetry group, N. So far, the statement of holography (or the 
gauge/gravity duality) is made about dualities between string theories and field the- 
ories in the limit of large N. On the other hand, we have used lattice simulations 
of SU(N) gauge theories with finite N (N=2,3). Therefore one may argue that the 
comparison we have made is not valid. However, the claim is that the discussed 
Wilson line does not depends strongly on N. This claim is based on results from 
recent lattice simulations of Wilson lines [1-3] for several values of N. It can also be 
seen in the lattice results we've been using (plotted in figure (3)), where there is no 
much difference between the SU(2) and SU(3) potentials. Assuming this claim to be 
correct, the large N Wilson line of the gauge theory should be similar to the ones we 
have used for the comparison. 

The second remark concerns quantum corrections. As the analysis in the present 
work is only at the classical level, the arising question is how quantum corrections 

^'^This requirement is motivated by the fact that the leading perturbative behavior of the QCD 
Wilson line is Coulombic. 

^^The KS potential at small distances is not exactly but rather approximately of the form ^ j. 
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would affect the results. In order to understand this point, we first review the 
situation in the flat space-time N.G. case, there the potential at the classical level is 
just a linear piece 



at all length scales. The full quantum expression (1.1) presents a dramatic change of 
this classical picture. While at long distances the quantum corrections are suppressed 
so the classical expression is a good approximation, at short length scales a tachyonic 
behavior takes over such that beyond some critical length (1.2) the full quantum 
expression (1.1) is not valid anymore. In the case of an asymptotically AdS curved 
background we expect the effect of the quantum corrections to be less dramatic, as 
we will explain below^^. In the classical level, the potential is linear at long distances 
but it becomes Coulombic (~ j) at short distances. The quantum corrections at 
long distances were proven to be suppressed; in [30] ([31], [32]) the effective action 
of closed (open) confining strings was found to obey the Nambu-Goto form up to 
high orders, i.e. deviations from the N.G. energy levels are suppressed as some high 
power of J (see also [33] where a discussion about quantum fluctuations of Wilson 
loop in string models is made). On dimensional grounds, the quantum corrections 
at short distances should have a Coulombic form (since the theory is asymptotically 
conformal) such that the potential would take the form 



The one loop correction to the classical term was calculated in [34] and the value of 
k was evaluated numerically in [35]^^ and analytically in [36] to be 



Therefore, the quantum corrections are suppressed in the limit where A — )■ oo so at 
short enough distances the classical expression will be dominant and no tachyonic 
behavior would pop out. We stress that this is true only at large A (or equivalently 
weakly curved background) which is not necessarily the case in QCD. To emphasis 
the difference between the fiat and curved space-time we remind that in the former 

^^However, there is an important difference between tire cases. While the tachyonic behavior 
in the case of flat space-time NG is an indication for an instability of the theory, in the present 
construction of the curved space-time case it is an indication for an instability of the specific Wilson 
line we chose and not of the theory. The reason for this is that in the latter case we have constrained 
the string's endpoints to end on the boundary while in the former there are still dynamical open 
strings. Therefore, by referring here to instability of the curved background case we only mean 
instability of the Wilson line and not of the theory. We leave this issue for future work. 

^^We thank Nadav Drukker for referring us to this paper. 
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case the expansion of the Wilson line is done in powers of ^ while in the latter the 
expansion is in powers of ^ (or equivalently in powers of a') due to the conformal 
symmetry. Hence, in the fiat space-time scenario the quantum corrections at some 
given a' would always become dominant at short enough distances and may present 
a tachyonic behavior, while in the asymptotically conformal curved space-time they 
are controllable. From this point of view we would prefer to consider the non-critical 
AdS^ background rather than the KS since the former is asymptotically AdS and 
the latter is not. 

Concerning the issue of quantum corrections, some future directions are pro- 
posed: 

• Following the statement made above, in asymptotically AdS backgrounds tachy- 
onic behavior can arise only at intermediate distances, as opposed to the fiat 
space-time case in which the tachyonic behavior is present below some critical 
length < £cr)- It would be interesting to characterize this finite range of 
tachyonic behavior at intermediate distances and in particular try to formulate 
conditions under which it is going to zero, such that no tachyons are involved. 

• Generalizing the quantum analysis to backgrounds which are not asymptoti- 
cally AdS (like the Dp-L^ branes). 
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A The conditions on the background to satisfy a coulombic 
potential in the UV 

This derivation is being done using notations and symbols of [5]. First we expand 
the /, g functions around u = oo 

f (u) = a^u^ + O [u^-^) 

g{u) = g (oo) + bju-^ + O {u-'^^) (A.l) 
k,j>0 

such that we can calculate the energy in the limit Uq ^ u\ using these expansions. 
The length of the string is given by 
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For the next steps we will need to use the invert function uq (i) 
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The fC function is given by 
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Using (A. 3) we can express it as a function of £ 
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Finally, the energy of the string as a function of its length would be given by 
E{i) = fiuo (£))-£-2/C {£) 
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(A.6) 

The conclusion is that in order to have a coulombic behaviour at small distances we 
have to demand the following asymptotic form of / 



ak 
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/ (m > Ma) ~ u^ 
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B Conical singularity in near extremal Dp brane 

We review the calculation of the periodicity of the compactified coordinate in the 
case of near extremal Dp brane background. The part of the metric which includes 
the radial coordiante, u, and the compactified one, x, is given by 



ds' 
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Following [37] let us define a new coordinate 

- = l+(^V (B.2) 

such that the singularity at m = ma maps to p = 0. Then for small p the relevant 2d 
part of the Euclidean metric (r = ix) is given by 

ds'^ ~ (7 - p) u^R-^p'^dr'^ + j^uf^R^dp^ (B.3) 

which can be written as [38] 

ds^ ~ n [dp^ + c^p'dT^] (B.4) 

where 

n = ^ R— (B.5) 

'^-P ii_p 2 



n 



In order to avoid a conical singularity the coordinate cr must be identified we a 
periodicity of 27r 

cr cr + 27r (B.7) 
means the periodicity of r is given by 

(3 = - = - R-^ B.8 

c 7 — p 

In a similar manner, the periodicity of the compactified coordinate of the non 
critical AdSn+i background was calculated in [39] to be 

n ua 

C Analytical expressions for energy and length of the string 
in near extremal Dp brane 

In the case of near extremal Dp brane, the metric is of the form 
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Inserting the above functions into the expression of the string's length (2.9) we have 
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we can write the length clS cL function of z 
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Notice that the length is measured in units of ua ■ ^ ^^^^ j ^ . Then we express it in 
terms of the universal system of units defined in (6) using the dictionary derived in 
section (7) (table (2)) to get 
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where dp = 3.4, 2.24 for p = 3,4, respectively. In a similar way we calculate the K, 
function (2.13) 
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Notice that the energy, which has the same dimension as /C, is measured in units of 
Expressing it in terms of the universal system of units we have 
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Finally, the energy is given by (2.12) 
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D The values of the constant terms in near extremal D„ brane 



In this appendix we will find the constant terms of the potential at short and long 
distances for the near extremal Dp brane, using the expressions for the energy and 
length derived in appendix (C) using the mass subtraction regularization scheme. 
Pay attention that we are before fixing the zero energy point as in [23] 's notations. 
First, at short distances, i.e. the limit z — > 0, one finds 
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where the three dots denote terms which vanish in this limit. The length is given by 
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such that 
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In fact, we could alreay conclude this power law behavior by the analysis in (A), but 
the important point here is the value of the constant term 
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Next we examine the behavior at long distances, i.e. z — )■ 1, there 

E = a-e-21C{z = l) 
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such that the constant term is given by 
Cl = -2/C {z = 1) 
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Finally, we conclude that the relative constant will be equal to 
AC = Cl- Cs 
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The expression in the square brackets appears to be equal to one such that 
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E Analytical expressions for energy and length of the string 
in non-critical AdSn+i 

Following the same steps as for the D branes example (C), the length of the string 
in the non-critical AdSn+i case is given by 
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Changing variables to y = z = (^^^ we have 
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Notice that in the present case length is measured in units of Then expressing 

it in terms of the universal system of units defined in (6) using the dictionary derived 
in section (7) (table (2)) we result with 
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In the same manner we derive the explicit expression for the /C function 
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We see that as in the D branes case, energy is measured in units of In the 
universal system of units the K, function will take the form 
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Finally the expression for the energy is given by 
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F The values of the constant terms in non-critical AdSn+i 

Following again the same steps as for the D brane cases in (D), we will find the 
constant terms of the energy in the cases of non-critical AdSn+i- We will use the 
expressions for the energy and length derived in appendix (E). Pay attention that we 
are before fixing the zero energy point as in [23] 's notations. First, at short distances, 
i.e. the limit 2 — )■ 0, one finds 
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where the three dots denote terms which vanish in this limit. The length is given by 
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We find that the constant term at short distances is given by 

a 5.O40Fr [1 - i] 
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In the opposite limit, 2 — )■ 1, corresponding to large distances, the energy takes the 
form 

E = a ■ i -2IC{z = 1) (F.6) 
such that the constant term is equal to 
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Finally, we conclude that the relative constant will be equal to 
AC = Cl- Cs 
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and this can be evaluate numerically for different values of n: 



AC = -5.34 



G The Wilson line of the Klebanov Strassler background 

The supergravity solution of the deformed conifold is of the following form [13, 27] 

ds'^ = h-^ {T)dxl^23 + hH^)d4 (G.l) 
where ds^ is the metric of the deformed conifold 
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K{t) = 

Therefore the / and g functions are 
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Near the boundary, i.e. at large r, the h (r) and K (r) function takes the form 
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Then the distance between the quarks is 
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Let us write this expression in the next form 
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such that the Kappa function will take the form 



/C(ro) 



-63' 



2i\/6 2^VQ 
Thus the energy in the region close to the boundary is given by 
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We would like to express the energy as a function of the distance between the 
quarks 
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To first order (at large tq) we have 
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such that the leading behavior of the energy is indeed Coulombic 
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It is not a surprise since what we have done by taking the leading behavior (G.19) is 
exactly considering only the leading geometry, which is AdS , as discussed in section 
(5). 

We now want to estimate the leading correction to the Coulombic potential. To 
do that we first find the leading correction to the dependence of the length on Tq by 
plugging in the leading order (G.19) into the subleaing term of (G.18) 
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Therefore the energy will take the form 
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Until now we have worked in the natural system of units of the KS model, but 
in order to make a comparison we have to switch to the universal system, as follows 
from the dictionary derived in section (7) (table (2)) 
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The constant term of the potential in the small distances regime can be read off 
from (G.23) 
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On the other hand, the constant term at large distances can be calculated numerically 
to be 



Cl = -2/C (r = 0) 
Therefore the relative constant is given by 
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